Josephson effect in superconducting wires supporting multiple Majorana edge states 
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We study superconducing-normabsuperconducting (SNS) Josephson junctions in ID topological 
superconductors which support more than one Majorana end mode. The variation of the energy 
spectrum with the superconducting phase is investigated by combining numerical diagonalizations 
of tight-binding models describing the SNS junction together with an analysis of appropriate low- 
energy effective Hamiltonians. We show that the 47r-periodicity characteristic of the fractional DC 
Josephson effect is preserved. Additionally, the ideal conductance of a NS junction with a topolog- 
ical supraconductor, hosting two Majorana modes at its ends, is doubled compared to the single 
Majorana case. Last, we illustrate how a non-zero superconducting phase gradient can potentially 
destroy the phases supporting multiple Majorana end states. 



I. INTRODUCTION 

Majorana fermions have recently attracted consider- 
able attention in the condensed matter community due 
to their exotic character (see recent reviews 1 and 1^. 
The fact that they obey non-Abelian statistics makes 
them interesting candidates for quantum computation.- 
Many systems hosting Majorana fermions have been pro- 
posed. One-dimensional semiconducting wires with a 
strong spin-orbit coupling in the proximity of a s-wave su- 
perconductor, and subject to a Zeeman field seem to offer 
a promising route to o bser ve Majorana fermions. Recent 
experiments using InStP^ or IuAJ^ quantum wires have 
reported signatures consistent with Majorana physics in 
conductance measurements in superconducting-normal 
(SN) junctions, though it has been shown that the in- 
terpretation of the data may not be so straightforward 
because of the presence of disorder-f^ltiSl 

Significant signatures of Majorana physics have also 
been predicted when two superconducting quantum wires 
are brought into contact to form a SNS Josephson junc- 
tion. It has first been shown by Kitaev^^ for single- 
channel topological superconducting wires that the two 
Majorana end states across the junction couple to each 
other and generate a fractional Josephson effect: instead 
of the usual 27r, the Josephson current exhibits a peri- 
odicity of Att with the phase difference between the two 
superconductors. This doubling of the periodicity can be 
interpreted as the tunneling of "half" of a Cooper pair, 
hence the term fractional. This prediction has later been 
extended to many different systems .1^^^^^ The fractional 
Josephson effect has been shown to be robust in realistic 
system, even when the Majorana fermions are coupled to 
other fermionic end states.l^Sl A direct consequence of the 
fractional Josephson effect is the presence of anomalous 
Shapiro steps at double the expected voltage; such dou- 
bling has recently been reported experimentally-I^S This 
fractional Josephson effect can be understood by a sim- 
ple low-energy effective Hamiltonian that couples the two 
Majorana fermions 71 and 72 across the ID junction, 
HoS = ir((/))7i72 where F is the coupling and the su- 



perconducting phase difference. The operator 17172 is 
obviously a conserved quantity of the Hamiltonian for all 
values of 0. This guaranties a level crossing at </) = tt in 
the evolution of the energy spectrum with the phase (f). 

However, the situation is far from obvious if we con- 
sider ID junctions between wires supporting several Ma- 
jorana fermions. The question we are addressing in this 
work is whether the fractional Josephson effect survives 
in such situations, and whether the periodicity of the 
phase/current dependence is modified. 

This question is relevant in the light of recent studies 
of a generalized ID Kitaev chain which supports two end 
Majorana fermions at each extremities, 'island of ID topo- 
logical superconducting wire that can support multiple 
end Majorana fermions. It has been recently shown by 
Tewari and SarP^that narrow semiconducting wires with 
a strong spin-orbit coupling, in the presence of a Zeeman 
splitting and in the proximity of an s-wave superconduc- 
tor belong to the symmetry class BDI in the classification 
of the non-interacting band Hamiltonians.^'* The BDI 
class implies a Z topological invariant instead of a Z2 
invariant for the D symmetry class.'23 Thus, a ID topo- 
logical superconducting wire belonging to the BDI class 
can support multiple end Majorana fermions. 

The main reason which can push ID non-interacting 
superconducting wires in the BDI class is the fact that 
in ID one can write a co mplete ly real Bogoliubov-de 
Gennes (BdG) Hamiltonian .122123] Along with the time- 
reversal symmetry (TRS) and an intrinsic particle-hole 
(p-h) symmetry of BdG Hamiltonians, there is a chi- 
ral symmetrjl2iH26| -^^hieh prevents the pairs of Majorana 
zero modes to recombine at the edge. The chiral sym- 
metry can be effectively realized in sufficiently narrow 
wires, typically for a wire width smaller than the spin- 
orbit coupling length.'2Slln such system, it was shown that 
two proximate robust Majorana fermions can be created 
by manipulating the sign of spin-orbit coupling and also 
that they are robust against disorder and small varia- 
tions of magnitude of magnetic fields-I^S! The presence of 
several Majorana fermions at one edge of the supercon- 
ducting wire opens several Andreev transport channels 
in SN junctions and therefore the conductance can reach 
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the value 2e2//i x Q with Q G 

The outHne of the paper is as follows: Sec. [lT| presents 
the tight-binding model for a chain supporting multiple 
Majorana at its ends, and some general symmetry ar- 
guments. Then Sec. Ill obtains the conductance of an 



NS junction where the superconductor is topological and 
supports two Majorana modes at each end. Section |IV| 
analyzes the energy levels and a low-energy effective 
Hamiltonian in Josephson junctions between topological 
superconductors with one and two Majorana modes at 
their ends: the 1 — 1 junction in |IV A[ the 2 — 2 junc- 
tion in [TVB] and the 1 - 2 junction in [TVCl Section |V] 
is devoted to breaking the Z phase through the addition 
of a uniform superconducting phase gradient. Section [Vl| 
presents the conclusions of the study. 



II. MODEL HAMILTONIAN AND SYMMETRY 
PROPERTIES 

In this section, we consider the ID model of spinless 
fermions, proposed by Niu et al. in Ref. |22j . which 
supports several Majorana end states. The Hamiltonian 
reads 



H 
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where Ai corresponds simultaneously to the nearest- 
neighbor (NN) hopping amplitude and to the nearest- 
neighbor superconducting gap, while A2 denotes the next- 
nearest-neighbor (NNN) hopping amplitude and next- 
nearest-neighbor superconducting gap. In what follows, 
the chemical potential /i is set to 1 and Ai is assumed 
positive. When A2 = 0, this Hamiltonian in Eq. ([I]) cor- 
responds to the Kitaev model. Before analyzing the 
properties of H , let us provide some general symmetry 
arguments which explain why a general ID Hamiltonian 
can sustain phases with more than one Majorana end 
states. 



Symmetry arguments and the role played by 
distant site hopping 



A ID superconducting system can have multiple Ma- 
jorana bound states at its ends when the system is 
time-reve rsal in variant (TRI) and exhibits particle-hole 
symmetryj ^^ l ^^ l For the two-band Bogoliubov-de Gennes 
(BdG) Hamiltonian presented here, this can be seen from 
the following simple argument. A general two-band BdG 
Hamiltonian H^dG obeys the particle-hole (p-h) symme- 
try by construction, and can be written in the p-h basis 
as 



BdG 



hr, 



(2) 



where r's are the Pauli matrices in the p-h space. Note 
that the p-h symmetry requires that 



hiik) = -h,{-k), 
h2{k) = -h^i-k), 
hsik) = hsi-k). 



(3) 



The time reversal operator for spinless fermions is just 
the complex conjugation operator. If the system is TRI, 
the components of h obey the following constraints 



hi{k)^hi{-k), 
h2{k) = ^h2{-k) 
hsik) = hsi-k). 



(4) 



Particle-hole and time-reversal symmetries impose the 
chiral symmetry represented by the operator ri which 
anti-commutes with the Hamiltonian 



{HBdG,Tl} = 0. 



(5) 



If -ffedG obeys these symmetries, then it follows that hi 
must vanish. Hence Hq^g has only two remaining com- 
ponents and h defines therefore a mapping from the Bril- 
louin Zone (BZ) to the Bloch "circle", h : ^ S\ The 
mapping is characterized by a winding number w which 
is an integer. Therefore a two-band BdG Hamiltonian 
belongs to the topolog ical BD I class characterized by a 
Z topological invariant .'2SI2II30] 

The expression of the winding number reads 



27r Jo 



dk{h X dkh.)i- 



The discrete form for the winding number reads 
w = -^ ^ sgn[h3dkh2]. 



(6) 



(7) 



Note from the symmetry constraints that the kernel of h2 
contains at least the special BZ points and tt. To cre- 
ate more Majorana bound states at one end, the winding 
number must have |?«| > 1. This implies that ker h2 must 
contain more points in addition to and tt. This can hap- 
pen by enlarging the unit cell through the addition of long 
range hopping. Then new nodes in the energy dispersion 
develop for k £ [0, tt] and can lead to higher winding num- 
ber. Hence the presence of higher-order hopping terms is 
a sufficient condition to have multiple Majorana bound 
states at the ends of a ID wire. 



B. Phase diagram 

The phase diagram of H has been established in 
Ref. [22] . Here we recover this phase diagram in a differ- 
ent manner, by using Eq. ([t]) to unambiguously charac- 
terize each topological phase in the (Ai,A2). This phase 
diagram is drawn for completeness in Fig. [l] 
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FIG. 1. (Color online.) Phase diagram for the Hamiltonian in 
Eq. ([l]). The insulating phases are separated by black lines. 
Each insulating phase is characterized by a winding number 
represented in red. 



The phase diagram associated with H is characterized 
by phases with w = 0, 1,2. Indeed, for A2 > 1 + Ai or 
A2 < — 1 and A2 < 1 — Ai, -ff can sustain a phase with 
two Majorana zero modes locahzed at each end.'^^l 



III. TRANSPORT IN SN JUNCTIONS 



where MiuS) = Tr[G^(w)rG'^(w)r(a;)], and G^(w) 
2[a;l + 2iV\^^ denotes the retarded Green's function. 
The differential conductance becomes 



dj[~uj + eV) df{uj - eV) 



which at T = reduces to 

dl _ 2e2 
dV ~ ~h 



duj 



M{eV). 



duj 



(9) 



(10) 



Taking an explicit trace over the transmission matrix, we 
find 



dl 8e2 



!det(r) 



dV h [(et/)2 -4det(r)]2 
which at zero bias becomes 

dl _ 4e2 



[2eyTr(r)]2' 



(11) 



(12) 



The zero-bias value of the differential conductance is 
thus the double of that expected for a junction with a 
single Majorana fermion at the interface. This is consis- 
tent with each interface Majorana contributing a jh 
to the total conductance. 



Before analyzing the physics of Josephson junctions 
made with wires supporting several Majorana fermions, 
we focus on an SN junction for which the SC is described 
by the Hamiltonian H in Eq. ([T]). For a junction be- 
tween a wire with one Majorana end state and a normal 
metal, it has been predicted that the differential conduc- 
tance exhibits a zero-bias peak of height jhW^ A 
similar question for a junction between a topological su- 
perconducting wire characterized by a topological index 
w > \ and a normal metal has been recently addressed 
by Diez et ai in Ref "2^]. The authors have shown that 
for such junctions the conductance G can reach a value 
of G = \w\ X 2e^ /h. Here we check this prediction by con- 
sidering a junction between a wire described by Eq. ([T]) 
supporting 4 Majorana fermions, 2 at each of its extrem- 
ities, and a normal wire. 

The low-energy properties of such a system are deter- 
mined by the coupling of the two end Majorana fermions 
with the normal metal. This coupling can be captured 
by a 2 X 2 hybridization matrix F. In order to com- 
pute the low-bias transport properties of this junction, 
one can directly use the S'-matrix formalism developed 
by Flensber^2ll^ noting that the two wave-functions 
for the Majorana fermions are orthogonal,'22l such that 
there is no inter-Majorana coupling term. 

The resulting expression for the current can be written 



IV. JOSEPHSON JUNCTIONS 

In this section, we analyze how the Jospehson effect is 
affected by the presence of several Majorana zero modes. 
Let us consider short Josephson junctions between two 
wires which can sustain several Majorana end states. The 
model Hamiltonian reads 

H = Hl+Hr + Ht (13) 

The Hamiltonian for the left wire, iJ/,, is described by the 
Hamiltonian H in Eq. ([T]) characterized by parameters 
(Af',A2). The right wire is characterized by Hn, which 
is obtained from Hl by changing (c|cj — >■ clcje*"^) in 
Eq. ([I]). Note that the same phase variable phase 4> is 
attached to the NN and NNN pairing terms. The second 
wire is characterized by the parameters (A{*,A|^). The 
tunneling Hamiltonian can be modeled as 

Ht = - ^(Afc]c,+i + Af c]_ic,+i -t- H.c.) (14) 
3 

Each wire is labeled by a topological index = 0, 1, 2 
with a — L, R. We have analyzed various — junc- 
tions which support Majorana fermions. We present be- 
low the 1 — 1, 1 — 2 and 2 — 2 junctions as sketched in 
Fig. [2] 



4 



7i 72 







71 



72 73 



7i 72 



73 74 



FIG. 2. To form Josephson junctions, wires characterized 
by winding numbers 1 — 1, 1 — 2 and 2 — 2 respectively, are 
brought into contact. Without loss of generality, the right- 
hand-side superconductors have real order parameters, while 
the left-hand-side superconductors have a SC phase <j). The 
low-energy Hamiltonian is assumed to contain only phase- 
dependent coupling terms between the Majorana fermions. 



Before analyzing the Josephson effect in junctions 
made with wires supporting multiple Majorana fermions, 
one may ask if a complex superconducting order param- 
eter in a topological superconducting wire can have an 
effect on its phase diagram. However, while Hii{(f)) is 
generically complex, an uniform phase (j) can be gauged 
away, yielding a real Hamiltonian and the same phase 
diagram as the one depicted in Fig. [T] 



A. The 1 — 1 Josephson junction 

Let us first consider a Josephson junction between 
two wires with a topological index it;" = 1. This 
type of junction has been extensively studied; here it 
is checked that the model yields results consistent with 
the known physics. The coupling constants are chosen 
to be (A",A2) — (1,1)- The numerical results for the 
dependence of the energy levels of this junction with the 
phase difference recover a 47r periodicity consistent with 
the anomalous Josephson effect (see Fig. [s]) . The eigen- 
value spectrum of the junction for (jj = 0,7r/2,7r is plot- 
ted in Fig. |3] Only for ip — it one recovers four zero- 
energy eigenvalues which corresponds to four Majorana 
fermions: one at each extremity and two at the junction. 

The 47r periodicity can be understood from a simple 
effective low energy Hamiltonian.'Sl The overlap between 
the wave functions of the two Majorana fermions at the 
extremities with the Majorana fermions at the interface is 
neglected. The si mplest low -energy effective Hamiltonian 
can be written sS^SMlEM 



H. 



off 



iti2 cos((^/2)7i72. 



(15) 



where 71 is a Majorana fermion at the right end of the 
first wire, 72 a Majorana fermion localized at the left end 
of the second wire and ii2 the effective tunneling ampli- 
tude (see Fig. [2]). One can check that the cosine behavior 
reproduces well the low energy spectrum. Introducing 
the fermion destruction operator c ~ ;^(7i + ^72)7 this 
is worth emphasizing that 17172 = c^c — ^ is trivially a 
conserved quantity of H^^^ . 
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FIG. 3. (Color online.) Top: The dependence of the lowest- 
energy eigenvalues on the SC phase difference between the two 
wires. Bottom: The eigenvalue spectrum for three different 
values of (j> — 0,7r/2,7r. The system exhibits four zero-energy 
modes at <^ = vr. The system parameters are Ai = 1 and 
A2 = 1. 



B. The 2 2 Josephson junction 

Let us consider now a junction between two topological 
superconductors characterized by winding numbers w" — 
2. 



1. Analysts of the spectrum 

The eigenvalue spectrum for a 2 — 2 junction is com- 
puted numerically. The result is shown in Fig. [4] The 
first important thing to note is that the anomalous 47r 
periodicity still holds. The only difference is that there 
are now four Majorana fermions forming at the junction 
when the phase difference is = tt. At each extrem- 
ity of the system there are also two Majorana fermions 
which subsist for any value of 0, making the ground state 
eight-fold degenerate at (/) — tt (see Fig. |4]). 

A low-energy Hamiltonian capable to describe the nu- 
merical results presented in Fig. |4] must involve four Ma- 
jorana fermions. Let us denote by 71,72 the Majorana 
fermions on the left side of the junction and by 73, 74 the 
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FIG. 4. (Color online.) Top: The dependence of the lowest- 
energy eigenvalues on the SC phase difference between the two 
wires. The energy branches are labeled using the occupation 
number representation for the two-fermion states formed from 
fusing four Majoranas across the junction. At low energy 
the two energy branches are linearly fitted: El'"^ ± E'^S^ = 
{d+iid-)\(j) — n\/2 and the value of coefficients d± is displayed. 
Bottom: The eigenvalue spectrum for three different values of 
4> = 0,7r/2,7r. The system exhibits eight zero-energy modes 
a,t (p = n. (Ai = 1, X2 — —1.5). 



Majorana fermions on the right side of the junction. A 
general Hamiltonian involving the four Majoranas can be 
written as 



2-2 

eff 



471(^1373 + ^1474) COS ■ 



suit of higher-order hopping processes between the two 
superconductors. Moreover, because they are odd func- 
tions of (j>, they are neglected in the following description 
of the low-energy spectrum at (f> ~ n. 

Contrary to the 1 — 1 junction, where there was an 
obvious quantity commuting with the Hamiltonian for 
all values of 4>, this quantity is not immediately evident 
in the 2 — 2 junction. In the later case, there are four 
Majorana fermions fusing to form two regular fermions 
that we denote by c± . This can be seen from the spectral 
decomposition of the effective Hamiltonian in Eq. (16) 



1) + 



1) 



where 



with d-i 



Ei 



d± cos(0/2). 



(17) 



(18) 



-- —K^ \/b^ ^JlP- — 4a^ with a — ti4i23 
2 V 2 

43 



ii3i 



13''24 



and h — t?4 + + t\^. Note that even if the 

/22 term in Eq. pG] ) was neglected, the form of the c± 
fermions as a function of the original Majorana fermions 
remains complicated. For example, under the reasonable 
assumption that the coupling strength between Majorana 
fermions situated at the same distance across the junction 
is identical, ^13 = ^24, if follows that 
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5± 
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1 



+ 24 
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: [74 -I- Z7i -I- g± (73 -I- 172)], with 



4i?3 + (tl4 



(19) 



(17 1 comutes triv- 



In this basis, the Hamiltonian in Eq. 
ially with the occupation numbers of the two fermions 
n± = Cj_c_|_. Therefore, the occupation numbers of the 
two fermions are conserved and one can use the eigenval- 
ues of the n± operators to label the energy branches of 
the two-particle states, namely the states as indi- 

cated in Fig. |4| . Such labeling of the four branches thus 
indicates that the level crossing is protected and hence 
guarantees the 47r-periodicity of the DC Josephson ef- 
fect for the 2 — 2 junction. Concerning the tight-binding 
Hamiltonian, the value of E'^^^ can be found numerically 
from a linear fit at low energy near </> = tt in Fig. [4j 



+i72(i2373 + ^2474) COS - 

+ («i347374 + «il27l72)/22 ('/'), 



(16) 



where /22 is an odd 27r periodic function of 0. The phase 
dependence of -ff^g ^ is fixed by enforcing a 27r periodicity 
together with the constraint that for = there is no 
direct coupling between 71 and 72, nor between 73 and 
74.12^ The cos(0/2) is required by gauge invariance. The 
effective Hamiltonian formally contains some direct tun- 
neling terms between the Majorana fermions on the same 
side of the junctions (^12 and ^34) which may be non-zero 
(/22(0) 7^ 0) when (/) ^ n-K. However, one can argue that 
such terms must be less significant since they are the re- 



2. Analysis of the Majorana polarization 

A useful tool to analyze the behavior of Majorana 
fermions is the Majorana polarization, a local topolog- 
ical order parameter introduced in Ref. [34j and related 
to the the degree of anomalous pairing in a ID topologi- 
cal wire. The existence of a Majorana fermion is recorded 
as a zero-energy Majorana polarization density of 0.5. A 
Majorana fermion can have a x- and j/-Majorana polar- 
ization. When the Hamiltonian is real, here at </> = and 
(j) = TT, the y-Majorana polarization Pm^ is 0. When the 
superconducting parameter acquires a phase, the polar- 
ization along the y direction is generally nonzero. 
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FIG. 5. (Color online.) Zero-energy Majorana polarization 
(a;-polarization in red and y- in blue) integrated over 40 sites 
around the junction. Four Majorana states form at the junc- 
tion for (j) = TT yielding a x-Majorana polarization of 2, and 
compensating the Majorana polarization of the four end MBS. 



Here, similar to Ref. [35] j the end Majorana fermions 
respond to the variation of (j> by rotating their Majo- 
rana polarization. This is due to the superconducting 
parameter becoming complex. However at = tt the 
superconducting parameter becomes real again and the 
zero-energy end states have the same Majorana polariza- 
tion. For this particular value of the phase, four new Ma- 
jorana fermions form at the junction, two in each wire, 
their x-Majorana polarization compensating the Majo- 
rana polarization of the end modes (see Fig. [5| . 



C. The 12 Josephson junction 

In the following, we restrict our attention to the 1 — 2 
junction. This is a particularly interesting problem since 
two distinct topological sectors are brought into contact 
via a Josephson junction. Here we assume that the left 
region is characterized by a winding number w — 1 and 
has the couplings Ai = I and X2 — 1, while the right 
region, with w — 2, has couplings Ai = 1, A2 = —1.5. 
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FIG. 6. (Color online.) Top: The dependence of the lowest- 
energy eigenvalues on the SC phase difference between the 
two wires. The dotted line indicate the fit of the lowest- 
energy nu mer ical eigenvalue with the analytical form in E+ 
from Eq. (21 1. The fit parameter is a = \/*i2 + ^la/^- The 
zero-energy state independent of (j) was not included in the 
figure. Bottom: The eigenvalue spectrum for three different 
values oi (j) = 0,7r/2,7r. The system exhibits six zero-energy 
modes at th — n. 



and 72,3 are proportional to cos(0/2) to compensate the 
sign change. A low-energy Hamiltonian describing the 
Majorana coupling in the 1 — 2 system can thus be written 
as 



1. Analysis of the spectrum 

It is expected that the physics of this junction to be 
dominated by three interacting Majorana fermions local- 
ized at the interface. One can show that for a winding 
number difference of one between the right and left side 
of the junction, one Majorana mode is bound at the inter- 
face for any choice of (f>. Moreover, at = tt the system 
is triply degenerate with three zero energy Majorana in 
the junction region (see Fig. |6]). 

Let 71 denote the Majorana fermion on the left side 
of the junction, and 72, 73 the Majorana fermions on 
the right side. The gauge invariance of the Hamiltonian 
suggests that the phase-dependent couplings between 71 



rrl-2 



«7l(^1272 



^1373) COS - -t- lt237273/l2(0), (20) 



where /12 is an odd 27r periodic function of 0. At = tt 
this Hamiltonian has three zero-energy eigenstates. For 
a phase difference of 7^ tt, the effective Hamiltonian in 



( 20 1 has one zero eigenvalue (required by the antisym- 
metry of the 3x3 matrix) and two non-zero eigenvalues. 
The constant zero-energy state, which is a Majorana edge 
state bound at the interface between two topologically 
nonequivalent regions, is the result of the difference of 
one unity between the topological indices of the two re- 
gions. 

Similar to the 2 — 2 junction, if the last term of Eq. pO| ) 
is neglected, the form for the two non-zero eigenvalues 
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FIG. 7. (Color online.) The Majorana polarization of the 
right-end Majorana fermions integrated over the last 30 sites. 
Note the oscillatory behavior of the x- (red) and y- (blue) 
Majorana polarization components with a total constant Ma- 
jorana density of 2 x 0.5, corresponding to two rotating Ma- 
jorana modes. 



becomes 



±2a cos((/)/2). 



tjg can 



The hopping dependent parameter a = j^/i^ 
be determined from a low-energy fit of the numerical dis- 
persion presented in Fig. |6] Note that taking into account 
also the term improves the quality of the fit, especially 
in the vicinity of the superconducting gap. Note also 
that, while the fit is accurate up to energies close to the 
superconducting gap, the effective Hamiltonian H^^^ is 
only a low-energy effective Hamiltonian. Therefore, one 
should not expect this Hamiltonian to describe the full 
dependence of the energy eigenvalues with the SC phase 
difference. 



FIG. 8. Zero-energy density of states integrated over the 40 
sites around the junction point. 
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FIG. 9. (Color online.) Zero-energy Majorana polarization 
integrated over the 40 sites around the junction. The x- 
polarization is depicted in red and the j/-polarization in blue. 
Three Majorana fermions having the same polarization are 
supported at the junction for (p = n. 



2. Majorana polarization analysis 

In the left-side wire, described by a winding number 
= 1, the superconducting parameter is chosen to be 
real and the single Majorana fermion at the left end is 
always fully x-polarized. For the right-hand-side wire, 
at the right end there are two Majorana fermions that 
respond to the twisting of the phase by gaining a y- 
polarization; The behavior of the x- and y-polarizations 
of these modes is presented in Fig. [T] 

Let us now analyze what happens at the junction be- 
tween the two wires. At (f> — tt, 3 Majorana fermions 
are expected. For any other value of 0, there is always 
at least one Majorana fermion stuck at the interface be- 
tween the two topologically nonequivalent regions. This 
can be seen by plotting the zero-energy density of states, 
as well as the zero-energy Majorana polarization, inte- 
grated over the 40 sites around the junction. In Fig. [8] 
we plot the integrated zero-energy density of states. Note 



that for (f> TT the density of states is constant and equal 
to 0.5 corresponding to a single Majorana mode bound 
at the junction. The sharp jump between 0.5 and 1.5 at 
(j) = TT describes the contribution of two extra Majorana 
modes which reach zero energy at 4> = n (these two extra 
zero-energy states appearing at (f) — tt can also be seen 
in the spectrum illustrated in Fig. [6|. This picture is 
confirmed by a plot of the zero-energy Majorana polar- 
ization, also integrated over the 40 sites around the junc- 
tion. The plot of the Majorana polarization is presented 
in Fig. |9] the jump in the a;-Majorana polarization at 
— TT can be understood as coming from the two Majo- 
rana modes that develop at zero energy. Thus, at = tt 
there are three Majorana fermions at the two extremi- 
ties of the wire (one at the left end and two at the right 
end), fully a;-polarized in the positive direction. These 
fermions are compensated by three Majorana fermions 
which form in the junction which are fully x-polarized in 
the opposite direction. 
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FIG. 10. Zero-energy local density of states. Addition of a 
small gradient V(/> = 0.01 per site leads to the removal of 
the protection of the zero-modes and a hybridization of the 
right-end Majorana fermions. 



V. WIRES WITH AN INHOMOGENEOUS 
SUPERCONDUCTING PHASE 

It was shown by Niu et al^^ that, for a system with 
two Majorana modes at each end, the four zero modes are 
present provided the time-reversal and chiral symmetries 
are not broken, namely as long as the system is placed 
in the BDI class .I^^IS^ Breaking time-reversal symmetry 
leads to the removal of the protection for the two addi- 
tional Majorana fermions since the system now belongs 
to the D symmetry class. Then the system can return 
to the more typical state with at most one Majorana 
fermion at each end. 

Here the fragility of the topological phase is explored 
again by breaking the TRS through the action of a uni- 
form superconducting phase gradient. Let us take a wire 
characterized by the topological index w — 2, with pa- 
rameters Ai = 1 and A2 = —1.5. A uniform phase gra- 
dient can be induced by the presence of supercurrents in 
the bulk of the superconductor. It was already shown 
that the phase gradient can be used to manipulate the 
creation or destruction of Majorana fermions-I^Sl If the 
entire wire is subjected to the phase gradient, the BdG 
Hamiltonian gains a TRS breaking /ii(fc) odd component: 

hi{k) = Ai sin(fc) sin(V0/2) + A2 sin(2A:) sin(V0), 

h2{k) = -Ai sin(fc) sin(V0/2) - A2 sin(2fc) cos(V(/)), (22) 

h^{k) = 2 - 2Ai cos(fc) cos(V(/)/2) - 2A2 cos(2fc) cos(V0), 

where the gradient V0 is the change of the phase over 
one site. Note that a phase gradient also creates a non- 
vanishing /iQ component that multiplies the identity Pauli 
matrix. Such a term tends to destroy the Cooper pairs. 
Let us assume a phase gradient, leading to an increase in 
the superconducting phase from left to right in a 100 site 
system. Because the Hamiltonian is now complex, the 
system is no longer characterized by a winding number 
w, but instead by the Kitaev Z2 invariant.'^ 

If the entire wire experiences a uniform phase gradient, 
then it is in a topologically nontrivial Z2 phase when the 
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FIG. 11. (Color online.) Zero-energy Majorana density inte- 
grated over 20 sites on the end left side (in blue) and on the 
right side (in red) of the wire. In the two regions predicted by 
by condition ( |23[ ) (which are delimited by the dashed lines for 
the chosen values of parameters Ai = 1, A2 = —1.5) there is 
exactly one Majorana fermion at each side indicating a non- 
trivial Z2 topological phase. However on the left side, there 
is a wide range of values of the gradient where there is still 
an unprotected Majorana doublet. 



condition 



|1 - A2 cos(V<^)| < |Ai cos(V(/)/2)| 



(23) 



is satisfied. Note that at vanishing phase gradient, one re- 
covers the phase diagram in Fig. [TJ with the only change 
that w — 2 and w — are both trivial Z2 phases. There- 
fore the addition of a small gradient Vcj) removes the 
protection of the Majorana doublets at the end of the 
wire. The accumulation of a SC phase near the left re- 
gion of the wire leads to a split of the zero energy modes 
into regular fermionic states (see Fig. 10 1. Note that for 
a strong enough phase gradient obeying the condition in 
the insulator makes a transition into a non- 



(23), 



Eq, 

trivial Z2 phase characterized by the presence of a single 
Majorana state at each end. This case is illustrated in 
Fig. |11| where we plot the zero-energy Majorana density 
integrated over the last 20 sites of the right and left side 
of the wire as a function of the gradient. The phase gradi- 
ent generally destroys the phase with multiple Majorana 
fermions since it renders complex the Hamiltonian which 
describes the wire. Note however that when the gradi- 
ent is close to or TT there is an unusual persistence of a 
doublet of Majorana fermions on the left side of the wire. 
This might be explained by observing that the supercon- 
ducting phase starts at zero on the left side of the wire 
and the superconducting parameters have only a small 
complex component for V0 ~ or tt. Then the Majo- 
rana doublet could be understood as a remnant of the 
Z phase for a small chiral symmetry breaking. Never- 
theless, the general conclusion remains valid: the phase 
gradient can serve as an external knob to switch from a 
Z topological superconductor to either a trivial or non- 
trivial Z2 insulator. 
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VI. CONCLUSIONS 

In this paper, we have studied in detail the Josephson 
effect between narrow semiconducting wires described by 
band Hamihonians in the BDI symmetry class which can 
support multiple Majorana end modes. Such situation 
may be encountered in relatively narrow semiconducting 
wires with a strong spin-orbit coupling, in the presence 
of a Zeeman splitting and in the proximity of an s-wave 
superconductor. 

Following Niu et al.^we mainly focused on two-band 
Hamiltonians resulting from nearest- and next-nearest- 
neighbor hopping integrals in a tight-binding description 
of such narrow wires. Based on exact diagonalization of a 
tight-binding model hosting multiple Majorana fermions 
together with a low-energy description of the junction, we 



have shown that the 47r anomalous periodicity of the DC 
Josephson effect is maintained. Even though the present 
work considered junctions between topological wires sup- 
porting one or two Majorana fermions at their extremi- 
ties, the low-energy description can be straightforwardly 
extended to junctions between wires supporting an arbi- 
trary number of Majorana fermions. 

Finally, we have shown that a uniform superconducting 
phase gradient is able to lift the protection for the states 
with multiple Majorana fermions and push the system 
into a non-trivial Z2 state characterized by the presence 
of a unique Majorana state at each end. 

The authors wish to thank J. Alicea and R. Aguado 
for stimulating discussions. The work of C.B. was sup- 
ported by the ERG Starting Independent Researcher 
Grant NANO-GRAPHENE 256965. 
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